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Thin sheets deposited on a substrate and interfaces of correlated materials offer a plethora of routes 
towards the realization of exotic phases of matter. In these systems, inversion symmetry is broken 
which strongly affects the properties of possible instabilities - in particular in the superconducting 
channel. By combining symmetry and energetic arguments, we derive general and experimentally ac¬ 
cessible selection rules for Cooper instabilities in noncentrosymmetric systems which yield necessary 
and sufficient conditions for spontaneous time-reversal-symmetry breaking at the superconducting 
transition and constrain the orientation of the triplet vector. We discuss in detail the implications 
for various different materials. For instance, we conclude that the pairing state in thin layers of 
Sr 2 Ru 04 must, as opposed to its bulk superconducting state, preserve time-reversal symmetry with 
its triplet vector being parallel to the plane of the system. All pairing states of this system allowed 
by the selection rules are predicted to display topological Majorana modes at dislocations or at 
the edge of the system. Applying our results to the LaAlOs/SrTiOs heterostructures, we find that 
while the condensates of the (001) and (110) oriented interfaces must be time-reversal symmetric, 
spontaneous time-reversal-symmetry breaking can only occur for the less studied (111) interface. 

We also discuss the consequences for thin layers of URu 2 Si 2 and UPt 3 as well as for single-layer 
FeSe. On a more general level, our considerations might serve as a design principle in the search for 
time-reversal-symmetry-breaking superconductivity in the absence of external magnetic fields. 


INTRODUCTION 

The realization and characterization of two- 
dimensional (2D) superconducting phases in various 
different systems constitutes a topic of great current 
interest [1-5] . This is motivated by the promising 
role played by 2D superconductors in the search for 
topological Majorana modes and related applications 
[6], by the gate tunability of the electronic properties 
[2, 4, 7], and by the fundamental interest in supercon¬ 
ducting transitions in reduced dimensions. Particularly 
interesting examples are given by LaAlOa/SrTiOa 
heterostructures, that show very rich electronic behavior 
[8], and single-layer FeSe on [001] SrTiC >3 with signifi¬ 
cantly enhanced transition temperatures compared to its 
bulk value [3] . This also motivates closer inspection of 
superconducting thin films of other correlated materials 
such as Sr 2 RuC >4 [9] and UPt 3 [10] which, in addition, 
promises to offer insights into the electronic properties 
of the bulk material. 

Noncentrosymmetric 2D superconductors form a par¬ 
ticularly important class since inversion symmetry is nat¬ 
urally broken in the practical realization of 2D systems: 
As shown in Fig. 1(a), both for an interface as well as 
for a thin layer on a substrate (B is vacuum) or in an 
asymmetric environment (B not vacuum, but A ^ B), 
inversion symmetry is broken; it can only be restored in 
the case of a thin layer in a symmetric environment as 
shown in Fig. 1(b). 

A pivotal property of superconducting states is their 
behavior under the inversion of the time direction. Not 
only does it determine the topological classification [11] 



Figure 1. From a symmetry point of view, 2D systems (yel¬ 
low) can be grouped into those realized in an asymmetric (a) 
or symmetric (b) environment. Layered materials consisting 
of weakly coupled sheets as shown in (c) are not included in 
our definition of 2D systems. 


but also essentially influences the electromagnetic and 
thermal response of these systems [12]. 

In this work, constraints on possible pairing states 
of noncentrosymmetric systems are derived that follow 
from the combination of symmetry and energetic argu¬ 
ments. These “symergetic” selection rules state that any 
superconducting order parameter transforming under a 
complex or multidimensional representation of the point 
group Q of the normal state necessarily breaks time- 
reversal symmetry (TRS). For 2D systems, it furthermore 
holds that this is only possible if Q contains a threefold ro¬ 
tation symmetry. Finally, if Q includes a twofold rotation 
symmetry perpendicular to the plane of the system, 
the component of the triplet vector along the axis of 
must vanish. These results hold under the assumption 
that (i) the energetic splitting E so of the Fermi surfaces 
is larger than the superconducting order parameter and 
(ii) that the superconducting phase does not break trans- 
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lation invariance. We emphasize that our notion of 2D 
does not include strongly anisotropic three-dimensional 
systems illustrated in Fig. 1(c). 

In the following, we will first present a proof of the 
selection rules stated above and then discuss the conse¬ 
quences for the time-reversal and topological properties 
of possible pairing states in several different 2D materials. 

RESULTS 

In order to decide which superconducting states are 
possible, we consider a system with pairing Hamiltonian 

4 h k^ k + l (j4 A fc (A) + H - c -) i 

(1) 

already taking into account assumption (ii). The 
fermionic creation and annihilation operators t/j k and tpk 
are TV-component spinors that describe the spin and or¬ 
bital degrees of freedom as well as potentially relevant 
subbands that result from the confinement along the di¬ 
rection perpendicular to the plane of the system. Cor¬ 
respondingly, the normal state Hamiltonian hk and the 
pairing function A*, in Eq. (1) are TV x TV matrices. Note 
that this approach and the following analysis goes be¬ 
yond the pseudospin description that is commonly used 
[13-15] for studying pairing in systems with spin-orbit 
coupling. 

To begin with the constraints resulting from symme¬ 
tries, let us investigate the transformation properties of 
A k under time-reversal and the elements g of the point 
group G of the normal state. The time-reversal operator 
is given by 0 = TIC with unitary T and AT denoting com¬ 
plex conjugation. Time-reversal acts on the pairing field 
according to 

A k A TA*_ k T T . (2) 

Under a point group operation g £ Q holds A& — g -x 
R^,(g)A R -i^ k R^(g), where R v (g) and Ri\g) transform 
vectors and spinors, respectively. 

To identify the order parameter, we expand the pairing 
field 

d n 

a=EEAA ( 3 ) 

n fi—1 

with respect to the basis of TV x TV matrix fields Xk/j 
transforming under the irreducible representation (IR) n 
of G • Here d n is the dimensionality of the IR n and g™ are 
complex-valued coefficients. Note that, before analyzing 
fluctuations, we first have to determine the form and, in 
particular, the symmetry properties, of possible order pa¬ 
rameters which is the central theme of the present work. 


TABLE I. Possible pairing states for a system with Ca v point 
group. It is indicated whether the phase preserves (y) or 
breaks (n) TRS. Here X and Y are continuous functions on 
the whole Brillouin zone transforming as the in-plane mo¬ 
menta k x and k y . For future reference, we also show the 
associated triplet component of the order parameter with oy 
denoting Pauli matrices in spin space. Although our analysis 
is more general, we here focus, for simplicity, on pairing states 
that transform trivially in orbital/subband space. 
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Including fluctuations will modify the behavior of physi¬ 
cal quantities in the vicinity of the phase transition. As 
will be seen below, the superconducting transitions we in¬ 
vestigate are always second order on the mean-field level 
and, hence, a Ginzburg-Landau (GL) expansion can be 
used to determine the candidate pairing states. Taking 
into account the usual orthogonality relations of IRs [16], 
the free energy T assumes the form 

d„ 

41= J [°]+EE tt "( T )l<l 2+J >4W ( 4 ) 

n [i=l 

with higher order terms £ O (M 4 )- The coefficient 
a no (T) that first changes sign determines the IR n = no 
of the order parameter. If the representation no is real, 
the TRS of the normal state implies that the matrix fields 
Xl 'T t can be chosen to be Hermitian (see Supplementary 

Information S.l) and from Eq. (2) follows 7 ?™ -A ±( 77 ^)* 
for 0 2 = =pl. As the global phase of the order parameter 
can always be absorbed by a 1/(1) transformation of the 
fields, we need at least a two-dimensional ( d no > 1) order 
parameter vector with a nontrivial relative phase to break 
TRS. As can be seen, e.g., in Table I, where all symmetry- 
allowed order parameters for the point group C^ v [16] 
are summarized, only the pairing state e(M) transform¬ 
ing under the 2D IR E breaks TRS. Note that this is 
different in the case of complex representations, where 
time-reversed partners transform according to different 
IRs. Consequently, we have to identify pairing states ei¬ 
ther in complex or in multi-dimensional IRs to obtain a 
TRS-breaking superconductor. 


Hmf = 
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Weak-pairing limit 

To deduce the consequences resulting from the en¬ 
ergetic assumption (i), it is convenient to diagonalize 
the free Hamiltonian h k by the unitary transformation 
%l>ki = Y a (4>ka)i fka that is made of its eigenfunctions 
4> ha satisfying h k cj) ka = £ ka 4> ka ■ Since h k is time-reversal 
symmetric, i.e. Qhf-Q^ 1 = h- k: we know that @4> ka is 
an eigenstate of h_ k with the same energy. The broken 
inversion symmetry at the interface along with spin-orbit 
coupling further imply that the Fermi surfaces are non¬ 
degenerate in the generic case. This implies for the wave 
functions that [17] 

<t>ka = e lv *0(f)_ k a, (5) 



Figure 2. Due to assumption (i), the inter-Fermi surface ma¬ 
trix elements (gray arrow) can be neglected as they couple 
states separated by energies of order E so . 


Microscopic derivation of the GL expansion 

Let us first focus on a single IR no with the associated 
interaction (g > 0) 


where the phase factors must satisfy the condition e llfk = 
as a consequence of 0 2 = q=l. 

The Hamiltonian can now be cast in the quadratic 
form H = | Y k using the Nambu spinor 

i’ka = (/fca’ilfca e-i ‘^) T - Tlie off-diagonal elements 
of the associated Bogoliubov-de Gennes (BdG) Hamilto¬ 
nian, characterizing the superconducting state, are given 
hy D ka b — (4>ka\A k T^ \(j) k b^. 

We now consider the weak-pairing limit (see Fig. 2) 
that implies that partners of a Cooper pair always orig¬ 
inate within a given Fermi sheet and not between states 
of different sheets. If this is the case, it holds 

Dkab = ^ka&ab- ( 6 ) 

Note, this assumption does not exclude frequently dis¬ 
cussed pairing states that are due to interband inter¬ 
actions. It merely requires that anomalous averages 
are made of the same quantum numbers as the normal 
state. In this weak-pairing limit, we immediately ob¬ 
tain the eigenvalues of the BdG Hamiltonian h^ dG as 
E ka = ±{e 2 ka + lA^I 2 ) 1 / 2 , i.e. |A fco | is the supercon¬ 
ducting gap on the Fermi surface. 

The behavior of A ka under point group operations fol¬ 
lows from inserting Eq. (3) and using that the wave func¬ 
tions of non-degenerate Fermi surfaces must transform 
as <j) ka = e lf>k R^{g)(f)R v (g)ka with phase factors e lf>k . We 
obtain that the basis functions := {(^kaWk^l^ka) 
transform under the same, a-independent, IRs as the ma¬ 
trix fields Xk/i - Thus, once we have found the IR no under 
which the pairing field A& transforms, along with the as¬ 
sociated order parameter vector r] n ° = ^r;"°, ■ • • , rpj° J , 
we also know the symmetry properties of the gap function 

A (7) 

11=1 

as it transforms exactly the same way. 


Hint = -g 






in the Cooper channel. One can write down the GL ex¬ 
pansion in the weak-pairing limit to all orders in the order 
parameter A ka formally expressed in terms of Fermi sur¬ 
face averages (|A fca | 2i ) a . As discussed in more detail in 
the Supplementary Information S.2, resummation shows 
that [??”] > 0 and, hence, the superconducting tran¬ 
sition must be second order on the mean-field level as long 
as the normal phase is time-reversal symmetric. This jus¬ 
tifies focusing on the first few orders of the GL expan¬ 
sion to deduce constraints on possible pairing states. To 
fourth order, it holds = I\r) n °\ A l3(z n °) + ..., where 
z n ° := Tj n °/\rj n °\, I > 0 is a (temperature-dependent) 
prefactor, and 


a \ /j,— 1 


(9) 


has been introduced. Let us assume that the minimum 
occurs at z n ° £ R d "o and define z?° := z"°| »« , ”o fi j 
for some po with 0. It follows from Eq. (9) that 


/3(z^°) ~ f3{z n °) — C(z n °)£ 2 (10) 


as £ —► 0. For no being a real and multidimensional IR, 
one finds C(z n °) > 0 (see Methods) except for z n ° = e Mo , 
with e M denoting the unit vector along the p direction, 
where C — 0 following from gauge invariance. In the 
latter case, one has to take instead z^° oc (e Mo + J 
with pi ^ po again yielding /3(z^°) < /3(z n °) for small, 
but finite £. This means that introducing relative com¬ 
plexity between the components lowers the free energy 
and, hence, the order parameter must break TRS. This 
means that the two pairing states eno) and in Ta¬ 
ble I are allowed by symmetry but suppressed energeti¬ 
cally in the weak-pairing limit. The analogous discussion 
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for complex IRs, which are best though of as real re¬ 
ducible representations of dimension 2 d nQ1 can be found 
in the Supplementary Information S.3. It yields that the 
superconducting state automatically breaks TRS when 
no is complex (also for d no = 1 ). 

Two remarks are in order. First, we emphasize that, al¬ 
though the main focus of this paper is on 2D systems, this 
result also holds for the three-dimensional (3D) case. Sec¬ 
ond, it provides a tool to identify TRS-breaking noncen- 
trosymmetric superconductors experimentally: The ob¬ 
servation of a splitting of the superconducting transition 
under the influence of an external symmetry-breaking 
perturbation such as uniaxial strain indicates that the 
representation of the order parameter must be multidi¬ 
mensional or complex implying broken TRS. 


Further consequences for 2D systems 

Shifting k —> — k in the weak-coupling limit of the pair¬ 
ing term and using the behavior of the phase factors in 
Eq. (5) under this shift, we obtain the important prop¬ 
erty 


A ka = ±A_ ka if 0 2 = =Fl- (11) 

Naturally, the upper sign is most relevant for fermionic 
pairing, yet we include the more general behavior for 
two reasons: Firstly, it illustrates the importance of 
normal state TRS for the fact that the gap function 
A ka = {(t) ka \A k T^\(j) ka ) has a well-defined parity. Sec¬ 
ondly, there are situations [18] where fermionic TRS is 
broken, however, the effective low-energy theory of the 
system has an emergent TRS that satisfies 0 2 = 1. 

Let us first focus on the upper signs in Eq. (11). Sup¬ 
pose that the point group contains a twofold rotation Gf 
with R v = — 1, i.e., perpendicular to the plane of the sys¬ 
tem. This is only allowed in even space dimensions, since 
det R v = 1, which is why we will be focusing on 2D sys¬ 
tems in the following. As dictated by Eq. (11), A ka has 
to be an even function of k and, hence, no solutions with 
finite gap can occur that are odd under this rotation. 

Before discussing in the next subsection under which 
conditions pairing states with vanishing intraband ma¬ 
trix elements are energetically disfavored, let us directly 
deduce the physical consequences of this selection rule. 
In the simplest case of just a single relevant orbital, the 
triplet vector d transforms as a vector under rotation 
forcing its component along the axis of G^~ to vanish 
in any pairing state that is even under C For the 
case of several orbitals/subbands, where the analysis is 
more involved (see Supplementary Information S.4), the 
orbital/subband-diagonal matrix elements of this compo¬ 
nent can be shown to vanish as well. 

Secondly, in order to discuss TRS, we start by con¬ 
sidering again the point group C± v as an example. The 


order parameter cannot transform under the 2D IR E 
that is required for TRS-breaking since E is odd under 
C 2 ■ Thus, we can exclude a finite order parameter that 
transforms as k x + ik y or any other superpositions of k x 
and k y for that matter. In the case of 2D systems, the 
matrix elements of the pairing field on a non-degenerate 
Fermi surface are too restrictive to allow for any of these 
pairing states. 

It is straightforward to generalize this analysis to all 
possible point groups of noncentrosymmetric 2D electron 
systems: For analogous reasons to G^., TRS-breaking 
superconductivity is not possible for the interface point 
group C4. The same holds for the isomorphic groups 
Di, D 2 d and S 4 that describe possible symmetries of 2D 
electronic sheets. For all other symmetry groups without 
any rotation symmetry (such as for Ci, i.e. in the absence 
of any symmetries) or containing only a twofold rotation 
normal to the plane, all IRs are real and one-dimensional 
such that TRS-breaking superconductivity is forbidden 
as well. For the remaining possible noncentrosymmetric 
point groups of 2D electron systems, all of which contain 
a threefold rotation, one cannot exclude TRS-breaking 
superconductivity without further assumptions. 

Similar reasoning implies that in the case of spinless 
fermions, i.e., for the lower signs in Eq. (11), the emergent 
TRS must be necessarily broken at the superconducting 
transition if the point group [19] contains a proper or 
improper fourfold rotation symmetry. 


Beyond weak pairing 


Let us next go beyond the weak-pairing limit and clar¬ 
ify under which conditions a (translation-invariant) su¬ 
perconducting phase with a vanishing intra-Fermi-surface 
order parameter can occur. 

To this end, we consider the most favorable scenario 
where, at low energies, the effective electron-electron in¬ 
teraction is dominated by the Cooper channel in Eq. ( 8 ) 
with no being odd under the twofold rotation G^~ £ Q. 
All other competing channels are assumed to be negligi¬ 
ble. From the analysis above, we know that the associ¬ 
ated matrix elements D ka i, vanish for a = b. By deriving 
a general upper bound for the energy gain A E when en¬ 
tering the superconducting phase at zero temperature, 
we have shown that superconductivity with D kaa = 0 
is only possible when the spin-orbit splitting E so on the 
Fermi surface satisfies 


E S o £ 


2A 


sinh(l/A) 


( 12 ) 


where A denotes the associated dimensionless coupling 
constant and A the cutoff of the theory (see Methods for 
more details). This means that, in the weak-coupling 
limit, A « 1, superconductivity can only emerge when 
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Figure 3. (a) Illustration of the gedankenexperiment explained in the main text with superconductor and substrate shown 
in yellow and green, respectively, (b) Fermi surfaces of thin layers of Sr 2 RuC >4 following from the model defined in the 
Methods, (c) and (d) show the low-energy part of the spectrum of the superconducting states transforming under A± and 
B 1 for open (periodic) boundary conditions along the y ( x ) direction. For concreteness, we have taken pure triplet pairing 
with d oc 0.2fi(sin fe, — sin ki, 0) T and d tx 0.08ti(sin fe, sin ki, 0) T , respectively (cf. Table I). As argued in the main text, the 
spectra of the B 2 and A 2 pairing states look qualitatively similar to (c) and (d). 


the spin-orbit coupling is exponentially small. The phys¬ 
ical reason is that the superconducting order parameter 
only couples states at energies differing by E so as illus¬ 
trated in Fig. 2. We have furthermore found that the 
value of the order parameter maximizing the energetic 
gain A E is larger than E so /2 which shows the relation 
between the validity of the weak-pairing description and 
assumption (i). 


Consequences for Sr^RuOi 

The general results presented above naturally lead to 
the gedankenexperiment illustrated in Fig. 3(a): Imagine 
putting a 3D bulk superconductor with a TRS-breaking 
order parameter on a substrate and gradually reducing its 
thickness d. If the resulting, necessarily noncentrosym- 
metric, point group Q of the thin layer system does not 
contain a threefold rotation symmetry, either supercon¬ 
ductivity disappears or TRS must be restored below a 
critical value of d. Furthermore, if C 2 £ Q , the triplet 
vector must be aligned parallel to the plane of the system. 
We note that these expectations are consistent with the 
explicit single-band-model calculations in Refs. 20 and 
21 . 

A natural candidate material for this gedankenexperi¬ 


ment is provided by Sr 2 Ru 04 hosting a superconducting 
phase [22] which is widely believed to be a TRS-breaking 
chiral p -wave state with triplet vector d & ex (k x + ik y )e z 
[23-25]. Its small superconducting transition tempera¬ 
ture and the near degeneracy [26] of the chiral p -wave 
order parameter with the triplet states transforming un¬ 
der the ID IRs of its bulk point group D 4 ^ along with 
the strong spin-orbit coupling of Ru make it an ideal 
system to apply the weak-pairing description and the se¬ 
lection rules derived above. Furthermore, thin layers of 
this material have been fabricated and shown to be su¬ 
perconducting [9[. As the point group is reduced to C± v 
by the presence of the substrate, we conclude that both 
TRS must be restored and the triplet vector must rotate 
to be aligned in the plane of the system upon reducing 
the thickness d. 

Consequently, the superconducting condensate belongs 
to symmetry class Dill [11] and is, hence, characterized 
by a 7Li topological invariant v with v = — 1 (y = + 1 ) 
characterizing the topological (trivial) phase [ 6 ] . In the 
weak-pairing limit, it assumes the simple form [27] 

v = II ( si g n (Afc a a)) , (13) 

a 

where the product involves all Fermi surfaces, k a is an 
arbitrary momentum on and m a the number of time- 
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reversal invariant momenta (green dots in Fig. 3(b)) en¬ 
closed by Fermi surface a. In order to calculate the 
topological properties of the remaining candidate pair¬ 
ing states transforming under the four ID IRs of C± v 
(see Table I) , we take the tight-binding model for the Ru 
t 2 g states with atomic spin-orbit coupling that is com¬ 
monly used [28] for bulk Sr 2 Ru 04 (see Methods) adding 
the inversion asymmetric hopping term S(L X sin(k y ) — 
L y sin(fc x )) allowed by the residual Cn, symmetry. Here 
Lj denote angular momentum operators projected onto 
the t 2 g manifold. The prefactor S is nonuniversal and 
unknown, however, the following discussion is indepen¬ 
dent of its value as long as all m a are the same as in the 
limit <5 = 0 which holds as long as S < 0.47ti with t\ de¬ 
noting the largest centrosymmetric hopping parameter. 
The resulting Fermi surfaces for the rather large value 
S = 0.45fi are shown in Fig. 3(b). 

In addition, we assume that the triplet component is 
larger than the admixed singlet component since triplet 
is dominant in bulk S^RuCH. Consequently, the sign of 
Afc aa is opposite on Fermi surfaces split by the inversion 
asymmetric hopping S ^ 0. From Eq. (13) and Fig. 3(b) 
it then directly follows that the state transforming un¬ 
der Ai, which is fully gapped in the weak-pairing limit, 
is topological. This is confirmed by the spectrum shown 
in Fig. 3(c) characterized by a gapless Kramers pair of 
Majorana modes crossing the bulk gap in the vicinity 
of k x = 7 r. The subgap states around k x = 0 result 
from the j3 and 7 band being topological separately. The 
band mixing, however, introduces a mass to the associ¬ 
ated edge modes. Since Eq. (13) also holds for the ID 
Dill Z 2 invariant (with m a = 1) [27], we directly find the 
nontrivial invariant u x = — 1 (v y = — 1 ) for the fictitious 
ID system at k x = 7r ( k y = 7 r). Therefore, the system 
is characterized by the weak indices ( 1 , 1 ) such that a 
Kramers doublet of isolated Majorana modes emerges at 
any dislocation with Burgers vector b = (&i, 62) satisfying 
h+b 2 odd [29, 30], 

The orbital polarization and symmetry restrictions 
render the Fermi-surface splitting very small along the 
high-symmetry lines k x = 0,7r and k y = 0,7r as can be 
seen in Fig. 3(b). Although being nodal in the strict 
weak-pairing limit, the B 2 pairing state is thus fully 
gapped for values of the order parameter that are larger 
than the small splitting at these high-symmetry lines but 
much smaller than the average value of E so such that the 
symergetic arguments presented above still apply. For 
this reason, the B 2 state has the same topological signa¬ 
tures as the A\ order parameter. 

However, the B\ and A 2 states have nodes along the 
T-M direction and, hence, the topological invariant v of 
the 2D system is ill defined. Nonetheless, the nontrivial 
fictitious ID invariant v x implies the presence of Majo¬ 
rana modes around k x = n at an interface parallel to 
the x axis which is confirmed by Fig. 3(d). This only 
holds as long as translation symmetry is preserved along 


x which might be irrelevant experimentally as very clean 
samples are already required to stabilize the supercon¬ 
ducting state itself. 

Taken together, this discussion implies that thin films 
of Sr 2 Ru 04 represent a promising system for the realiza¬ 
tion of Majorana modes. 


Other TRS-breaking condensates 

Let us discuss two further materials, URu 2 Si 2 and 
UPt 3 , which are believed to host a TRS-breaking super¬ 
conducting bulk state [31, 32]. 

To begin with URii 2 Si 2 , we first note that it is from a 
symergetic point of view very similar to S^RuCH: The 
point group D^h of its bulk is reduced to C^ v in a thin film 
of (001) orientation and the combination of the strong 
spin-orbit coupling of U and the small transition temper¬ 
ature [33] make the weak-pairing description possible. It 
follows that, if the thin film still displays superconduc¬ 
tivity, TRS must be restored in the condensate and the 
triplet component of the order parameter must be aligned 
parallel to the plane of the system. In particular, from a 
pure symmetry point of view, the most natural [34] can¬ 
didate pairing state transforming under the IR E 
subduced from E g of D^h, which is the IR of the bulk 
pairing state A& = ia y (k x +ik y )k z [31, 35], is suppressed 
in the weak-pairing limit. Consequently, superconductiv¬ 
ity is likely to disappear in the thin-layer limit. 

As for UPt 3 , the weak-pairing approximation is ex¬ 
pected to be applicable for the same reason as in the case 
of URu 2 Si 2 and S^RuCH, but the point symmetries are 
different: The bulk point group D 6 h is reduced to Cq v in a 
(001) film and, hence, contains a threefold rotation sym¬ 
metry such that TRS-breaking cannot be excluded. Due 
to C% € Cq v , the triplet vector which is largely aligned 
along the z axis in the bulk condensate [36] must rotate 
to be parallel to the xy plane. Note that superconductiv¬ 
ity has been reported in (OOl)-oriented films of UPts [10]. 
From a symmetry point of view, there are 10 possible su¬ 
perconducting phases - four associated with the four ID 
IRs of Cq v and 3 with each of the two 2D IRs E\ and E 2 . 
From the symergetic arguments presented above we ex¬ 
clude Ei and two of the ID IRs as these representations 
are odd under C 2 - Furthermore, the two time-reversal 
symmetric order parameter configurations transforming 
under E 2 can be discarded and we are left with only 3 
candidate pairing states for the thin layer system: The 
time-reversal symmetric s-wave and i-wave states trans¬ 
forming under A\ and A 2 , respectively, as well as the 
TRS-breaking state transforming as (X + iY) 2 . 
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Further examples 

Finally, let us discuss two additional 2D superconduct¬ 
ing systems without TRS-breaking 3D analogue. We be¬ 
gin with the LaAlOa/SrTiOa heterostructures that show 
interface conductivity [37-39] for the three different ori¬ 
entations (001), (HO), and (111) of the interface with re¬ 
spective point groups C± v , Ci v , and C^ v while supercon¬ 
ductivity has so far only been reported for the former two 
orientations [1, 5]. Due to the strong spin-orbit splitting 
of the Fermi surfaces [40, 41], the weak-pairing descrip¬ 
tion is clearly appropriate. From our general symergetic 
arguments, it follows that the condensates of the (001) 
and (110) interfaces must be necessarily time-reversal 
symmetric whereas the (111) heterostructure allows for 
exotic TRS-breaking superconductivity. Due to the ab¬ 
sence of a C 2 symmetry it is also the first system we have 
discussed so far that makes an out-of-plane triplet vector 
possible. Taken together this motivates a closer exper¬ 
imental inspection of the low-temperature properties of 
the (111) interface. 

In order to calculate the topological invariant v in 
Eq. (13) of the (001) and (110) interfaces, a microscopic 
calculation has to be performed since there is no 3D ana¬ 
logue to compare with and the symmetry properties alone 
do not determine v. The analysis of Refs. [42, 43] shows 
that the topological properties are directly related to the 
origin (electron-phonon/purely electronic) of the interac¬ 
tion driving the superconducting instability. 

Our final example is single-layer FeSe on SrTiOs. If 
the weak-pairing description is also valid for this system, 
the symergetic restrictions apply and the superconductor 
cannot transform under the 2D IR of C^ v , thus, preserv¬ 
ing TRS. In combination with experiment indicating the 
absence of nodes [44] , there are only three possible pair¬ 
ing states: The pairing field can have the same sign on 
all four (spin-split) electron pockets around the M point 
(s+ +++ ), the signs can be pairwise identical (s +H ) 
or only differ on one Fermi surface (s ++H ). Only the 

latter pairing state has a nontrivial Dill invariant v as 
readily follows from Eq. (13). It has recently been shown 
[43] under very general assumptions that an s -1 " -1 " -1 state 
is not possible irrespective of whether superconductivity 
arises from the coupling to collective particle-hole modes 
or from phonons. Therefore, FeSe is most likely a topo¬ 
logically trivial, TRS-preserving superconductor. 


DISCUSSION 

The applications of the symergetic selection rules to 
various materials discussed above show that these can 
both be used to pinpoint the order parameter of 2D su¬ 
perconductors as well as serve as design principles in the 
search for superconducting phases with exotic properties 
such as broken TRS or nontrivial topologies. In this con¬ 


text, it is particularly important that our arguments are 
only based on the assumptions (i) and (ii) and, hence, go 
beyond model studies, i.e., do not depend on microscopic 
details such as number and character of relevant orbitals 
or the form of the interaction driving the superconduct¬ 
ing instability. 

Since inversion symmetry is locally broken at the sur¬ 
face of a material, one might wonder whether the symer¬ 
getic selection rules are also relevant for the behavior 
of the superconducting phase at the boundary of the 
system. In the case of a material such as Sr 2 Ru 04 
which consists of weakly coupled layers as illustrated in 
Fig. 1(c), the condensate near a surface perpendicular 
to the 2 axis can be thought of as a set quasi-2D sys¬ 
tems with E so increasing as the distance to the surface is 
reduced and, hence, bears strong similarities to the su¬ 
perconductor in our gedankenexperiment. In combina¬ 
tion with the near degeneracy [26] with the triplet states 
transforming under the ID IRs, we expect the triplet vec¬ 
tor to rotate to be parallel to surface and TRS to be re¬ 
stored locally. It is an interesting open question whether 
this could account for the absence of magnetic signals in 
experiment [45] that are expected from the chiral p -wave 
nature of the bulk order parameter. 


METHODS 

Fourth order of the GL expansion 

Let us provide more details on the proof by contradiction based 
on the fourth order GL expansion first focusing on real IRs. Ex¬ 
panding Eq. (9) with z n ° —> = z n ° | n 0 n 0 i( to leading 

t z v-o~' z Mo e 

nontrivial order in one finds Eq. (10) with 

C(Z-O) = 4 £ (« 0 ) 2 Y,{«T°< n a) 2 ) a ■ ( M1 ) 

a 

To derive Eq. (Ml), it has been taken into account that 0 £ R 
and that the symmetries require the free energy to be invariant 
under 

V/zo € {1,2,...,d no }, (M2) 

for any real IR no of 2D and 3D point groups. From Eq. (Ml), we 
directly see that C(z n °) > 0 for all z n ° / e M0 unless (p^ n o = 0. 
In the latter case, however, the superconducting state is fully un¬ 
gapped in the weak-pairing limit and, hence, disfavored energeti¬ 
cally as discussed in the main text. 

If z n ° = e M0 , we will use z^ = (e Mo + i^e^/y/l + £ 2 again 
yielding Eq. (10), but with modified 

C(z-O) = [(wr°)l - (WVD 3 )>J . (M3) 

a 

which is readily shown to be positive (as long as are not 

identically zero). This completes the proof for the case of real IRs 
of point groups. 

Due to TRS, complex IRs are always degenerate with their con¬ 
jugate IR and, hence, can be seen as reducible representations of 
doubled dimension. Being reducible, symmetries are less restrictive 
in this case and, in particular, Eq. (M2) is not guaranteed any more 
which necessitates a generalized form of the proof presented above. 
The latter can be found in the Supplementary Information S.3. 


f -»?M 0 , P = MO 



Inequality for the condensation energy 

Here we discuss how Eq. (12) of the main text is obtained. To 
derive a necessary condition for the emergence of a superconducting 
state with Dk aa = 0, we analyze whether its zero-temperature 
condensation energy 

1 N A 2 

AE(Ao) = - E E (l- E fca(Ao)l ^ kfcal) - ^ (M4) 

Z a=1 k Zg 

is positive for some finite Ao = g|? 7 no |. Here e*. a and E^ a denote 
the different bands of the normal state and of the superconducting 
mean-field Hamiltonian, respectively. 

To focus on the essential part of the physics, let us consider only 
N = 2 singly degenerate bands. Replacing |-Dfci 2 | by its maximum 
value Ao m yields an upper bound AE max (Ao) on the condensation 
energy. Physically, it corresponds to the situation of “optimal basis 
functions” with \D^\ 2 \ being constant except for negligibly small 
regions where it has to vanish as dictated by symmetry. Evaluating 
the sum in Eq. (M4) as an integral (cut off energetically at A, 
constant density of states pp) shows that the condensation energy 
can only be positive when the spin-orbit splitting E so on the Fermi 
surface satisfies Eq. (12), where the dimensionless coupling constant 
is defined by A = 2 ppm 2 g. Furthermore, one finds that Ao m > 
E so /2 at the positive maximum of AE max (Ao) revealing the direct 
connection to assumption (i). 


Model for Sr 2 Ru 04 


To be self contained, we define the model used in the main text 
to calculate the spectrum of Sr 2 Ru 04 . The centrosymmetric part 


hi = 


€xy (fc) P> &€xy 
0 
0 

\ E L r 

j=*,y, z 


0 


0 


£iz(fc) — y t ,, sin(fci) sin(fc 2 ) 


trj sin(fci) sin(fc 2 ) 


z(k)~ fi 


(M5) 


is taken to be of the form usually applied (see, e.g., 
Ref. 28) to describe the bulk of the material. In Eq. (M5), 
the orbital basis {4d xa , 4d xz , 4d yz } of Ru orbitals is used 
and <?j are Pauli matrices representing spin. Furthermore, 
e xl ,(fc) = -2f 3 (cos(fci) + cos(fc 2 )) - 4f 4 cos(fci) cos(fc 2 ), e xz (k) = 
— 2fi cos(fci) — 2t 2 cos(fc 2 ), and e yz (k) = —21 2 cos(fci) — 2ti cos(fc 2 ). 
Adding the inversion antisymmetric hopping term = 

S(L X sin (k y ) — L y sin (A,:)) already introduced in the main text de¬ 
fines the normal state Hamiltonian h k = h k + h k in Eq. (1). To 
obtain Fig. 3(b-d), we have taken t 2 = O.lti, t 3 = 0.8ti, t 4 = 0.3ti, 
ir; = —0.04ti, A = 0.2fi, fi = 1 1 and 5e xy = O.lti as deduced in 
Ref. 28. 
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S.l Hermiticity of basis functions 


Here we show that the basis functions {x] 1 , } introduced in Eq. (3) to expresses the superconducting order parameter 
can always be chosen to satisfy the Hermiticity constraint 

(xl,T^= X l^ (SI) 

as a consequence of the TRS of the normal phase. Here fi denotes the complex conjugate representation of the IR 
n. Eq. (SI) is central for the analysis of the main text as it determines the transformation behavior of the different 
superconducting states under time-reversal. It also leads to the property (^")* = Pka the weak-pairing basis 
functions The following proof of Eq. (SI) generalizes the results of Ref. 46 to complex representations and to 

multiband systems, i.e., goes beyond the pseudospin picture. 

We first generalize the parameterization (3) to 


Ak 


EE»; 


n^n 
k/i 


(K- 


K)t , 


+ hi 




Vu,hi,h% G <D. 


(S2) 




It is straightforwardly shown that 5^ ( (/i". h 2 ) transforms exactly as Xk^T* under Q by taking into account that 
[0, Rp(g)] = 0 for all g G Q. Upon reparameterization H 2 = 2 (/i" + h 2 ) and H 2 = 2i (h] 1 ' — h 2 ), it holds 


= E(H2)*). 

Applying Eq. (2) to Eq. (S2), one finds that time-reversal is represented by 

A (±(<)*,(A)*,(Ar) 

for 0 2 = =pl. The free-energy expansion has again the form (4) with 

2 

a n {T-H?,H2)= ]T Mj™) (T)(Hj‘)*H™, 

j,j '=i 


(S3) 


(S4) 


(S5) 


as follows from gauge invariance and basic scaling arguments. Reality of the free energy and invariance under time- 
reversal (S4) forces M^ to be Hermitian and = (Afy")) r , respectively. It implies that ( H n )* = H n (without 
loss of generality) at the minimum of the free energy and, together with Eq. (S3), yields Eq. (SI). In addition, we 
find from Eq. (S4) that time-reversal simply becomes 


for 0 2 = =Fl a s used in the main text. 


% 



(S6) 


S.2 Microscopic GL expansion and resummation 


The GL expansion is straightforwardly derived by decoupling the interaction in Eq. (8) via a Hubbard-Stratonovich 
transformation in the Cooper channel, integrating out the fermionic degrees of freedom in the weak-pairing limit and 
subsequent expansion in the superconducting order parameter. One finds 


where 


1 “"0 00 o2Z-l , 0 2Z' 

m°\ = m + - E K °\ 2 +E —(-i)^(a ,t) ^ (| E 

fi=l 1=1 a \ fi= 1 


h(A,T):=p F d 


-A 


A + A 


> o 


(S7) 


(S8) 
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with ui n representing fermionic Matsubara frequencies and (...) a denoting the average over Fermi surface a defined 

by 

(.. ,) a := pp 1 f dfl k Pai^lk) ■■■ ■ (S9) 

J a 

Here pp and p a (H) are the total and angle-/Fermi-surface-resolved density of states, respectively. Furthermore, f d flk 
in Eq. (S9) describes integration over Fermi surface a. While the explicit form of p a (H) is irrelevant for the derivation 
of the selection rules, we only take advantage of positivity /? a (H) > 0 . 

To study the order of the mean-field transition, we have to focus on the fourth and higher order terms J -> 4 of the 
GL expansion which we write as 



ri , 0 un.0 
Ip H'ka 



(S10) 


Evaluating Ii in Eq. (S 8 ) (for simplicity in the limit A — > oo) leads to a power series representation of the dimensionless 
function /> 4 (:r) that is positive but only converges for |a;| < 7 r 2 / 4 . 

To access larger values of the superconducting order parameter, one can interchange the summation over l and the 
integration with respect to energy e in Eq. (S7) yielding the analytic continuation on R + of the aforementioned series 
representation 

/> 4 (z) = T ~ 1 J ; g>i{y) = * (4y- ln(l + 4y)) . (Sll) 


Due to g >4 > 0 and, hence, /> 4 > 0 on R + , it follows -F> 4 [ 3 ™ 0 ] > 0 as stated in the main text. 


S.3 Complex IRs 

As mentioned in the main text, complex IRs n ^ n deserve further investigation as far as the proof by contradiction 
based on the fourth order GL expansion is concerned. Due to the representation (S 6 ) of time-reversal, the TRS of 
the high-temperature phase implies a n = a fi in the GL expansion (4). Consequently, n and h are always degenerate 
at the quadratic level and are more conveniently treated as a real reducible representation n of dimension 2 d n with 
basis functions 

Xlp = Ukp + Xlp) , Xlp +dn = {xlp, - Xkp) ( S12 ) 

and associated = {(^kahk^^ka) G H. 

Using the real representation associated with a pair of complex conjugate representations, one can repeat the 
proof by contradiction presented in the main text. However, this does not work for all pairs of complex IRs of the 
point groups of crystalline 2D and 3D systems: For the complex IRs of C 4 , 64 and , all of which are ID, the 
associated 2D real representation does not ensure invariance under Eq. (M2). Here, the symmetries only impose the 
constraints ((^feo) 4 ) a = a independent of p, = 1,2 and 712 = -721 = 7 with 7 ^/ = J2 a ((Vka) 3( Pka)a In serting 
= (cos((/>)e l £, sin(0)) in Eq. (9) and expanding in £ yields Eq. (10) with C = /3/i( 0,7//3) where 

fi{(f>,gi) = sin(20) (sin(2</>) + 31 cos(2(?!>)) (S13) 

and /3 = ]U) a ((^feaV’Ia) 2 ) a - F°r an y nonzero 31 , one can find <f> e R with /i(<(>, 31 ) < 0 indicating that the proof by 
contradiction of the main text has to be extended in order to exclude real-valued order parameter vectors. 

To this end, we consider a generalized transformation z£ —> e l ^' T z£ with r = {t x ,t v ,t z ) denoting Pauli matrices. 
So far, we have focused on £ = ^ 3 . Choosing £ = £e 1 instead, we obtain C = ^ 7 ^ 27 /( 0 ; — /?)) where 

/ 2 (</>, 32 ) = cos(2 (j>) (cos(2 cj)) + 32 sin(2^>)). (S14) 

One can show there is no (f> G R with both /i(<(>, 31 ), being non-positive unless 

gigi > 1 ^ 7 2 >^P(a- /?)• 


(S15) 
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Eq. (S15) can be shown to be not satisfied by noting that ) := ]T) a {fkaPka) defines an inner product and, 

hence, obeys the Cauchy-Schwarz inequality. The latter leads to 


7 2 = 


<^V S ] (V 1S ) 2 - (V**) 2 ) 


2 < \ </V"l^V”) • <(^") 2 - (7> 2 ") 2 |(/ K ) 2 - (<^) 2 > = \p{<* - /?). (S16) 


Consequently, in the case of a complex IR, the superconducting order parameter must necessarily break TRS if the 
weak-pairing description applies. 


S.4 Orientation of the triplet vector 

Finally, we present more details on the symergetic restrictions on the triplet vector in the case of several orbitals 
and/or subbands. In general, the decomposition into singlet and triplet assumes the form 

A k = A s k T+ Y, D i°i T ( S17 ) 

j=x,y,z 

with both A k and D k being matrices in orbital/subband space. For concreteness, we choose a real orbital basis for 
which T = 1 (g> itjy. Fermi statistics implies 


(A_fc) T = A k (dQ 


= ~Dl 


(S18) 


Focusing on the relevant situation C'/~ € G, the order parameter must be even under C% and, hence, satisfy 


R^)A b : k R\ o {G^) = A», R^ 0 {C^)D^ k R\ o (.C^) = -D^, R^{Ct)D- k R\ o {C^) = D A k (S19) 


R > 2 


1,2 




with Ry, 0 (C'/ _ ) denoting the representation of C 1 / - in orbital space. 

To proceed further, we take advantage of the fact that the eigenspaces of the spin-independent part of the unit cell 
Hamiltonian are spanned by basis functions of the IRs of Q (including spin-orbit coupling would require considering 
the IRs of the associated double group). For example, in case of S^RuCH the relevant low-energy subspaces are 
spanned by the 4 d xy ( B 2 of C 4 V ) and 4 d xz , M yz (E of C 4 V ) orbitals. This basis is most convenient for discussing the 
consequences for the triplet component of the order parameter as Rip 0 {C^~) simply becomes 

Ri’o( C 2 ) = dia S (ci,c 2) ... ,c N/2 ) , Cj <={+1,-1}, (S20) 


in this basis. Using this in Eq. (S19) in combination with the Fermi constraint (S18), it directly follows that all 
diagonal components of D k vanish identically as stated in the main text. 




